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Abstract 

We find conditions required to achieve complete population transfer, via coherent population 
trapping, from an initial state to a designated final state at a designated time in a degenerate 3- 
level atom, where transitions are caused by an external interaction. Complete population transfer 
from an initially occupied state 1 to a designated state 2 occurs under two conditions. First, there 
is a constraint on the ratios of the transition matrix elements of the external interaction. Second, 
there is a constraint on the action integral over the interaction, or "area", corresponding to the 
phase shift induced by the external interaction. Both conditions may be expressed in terms of 
simple odd integers. 
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I. INTRODUCTION 



Population control in quantum systems, namely transfer of electrons from an ensemble 
of atoms all in the same initial state to specified final states, is used in problems ranging 
from coherent population trapping 



transparency 6 L 3, 



dynamics 



B3, 



inc 

n 



and quantum computing [10, 



uding electromagnetically induced 



11 



12 



13 



14. 



3Q, 



to chemical 



17|, [l8|. These problems are modeled in terms of an n-level (often 3-level) atom 



interacting with a strong external field 



19 



20 



21 



22 



23 



24j . In this paper we consider 



population control in a nearly degenerate 3-level atom. We show that the 3-level atom 
is relatively easy to understand in the degenerate limit, where all 3 states have the same 
energy. We specifically show how to achieve complete population transfer in a degenerate 
3-level atom where the matrix elements of the external interaction have a common time 
dependence. 

Our degenerate a ppr oximation j^j is mathematically similar to the rotating wave ap- 
proximation (RWA) which has been widely applied to both 2 and 3-level atomic 
models. However, in the RWA degenerate atomic states are not used. Instead one tunes the 
frequency of the external field to the frequency difference of two non-degenerate levels so 
that the detuning parameter tends to zero. Thus, in RWA the initial state of an atom 
plus one photon is degenerate in energy with the final state of the atom. One advantage of 
using degenerate atomic states, as done in this paper, is that one may use external interac- 
tions with a broad range of frequencies. Another advantage is that the interaction frequency 
can be used as a control parameter, e.g. to vary the duration of time that the transferred 
population remains in the designated state, or to reduce the population leakage that occurs 
when the three levels are not fully degenerate. 

In the next section we derive analytic formulae for the probabilities, Pi^^it), that the 
electron is in level 1, 2 or 3 at time t. Then we seek conditions for complete population 
transfer to a designated level that is initially unoccupied at a designated time to. This places 
conditions on both the relative strengths of the interaction matrix elements, Vij(t), and the 
action integral, A(t ) = f Q ta V{t')dt' . We show that these conditions may be expressed in 
terms of two odd integers, ri\ and rii- Some analysis is done for population leakage, when 
the levels are not quite degenerate. Calculations are presented. 
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II. THEORY 



Let us consider an n-level atom interacting with an external field, V ex t(r, t). The total 
Hamiltonian for this system is H = H + V ext (f, t). The n eigenstates, <p k , and corresponding 
eigenenergies, E k , of H are assumed to be known. The total wavefunction may be expanded 
in terms of the known eigenstates, namely, ^{t) = ai(t)0i + ct2(^)02 + • • • + a n (t)<f) n . With 
atomic units, using = (H + V ext (r,t))^ , with H (f) k = E k <p k and J 4>*(f) k df = 6j k , one 
then obtains [soj ] . 

n 

idj(t) = E^t) + V jk (t)a k (t) , (1) 

k=l 

where Vj k (t) = J c()*V ex t{r,t)(f) k dr. These equations are exact for an n-level atom. 

We now require that the system be degenerate, namely that all the energies, Ej, are the 
same. Since the zero point of energy is arbitrary, one may generally set Ej = 0. These 
conditions give the coupled equations for a degenerate n-level system, namely, 

n 

iaj{t) = ^V jk {t)a k {t) . (2) 

k=l 

We use the initial condition a^O) = 1, and aj(0) = for j ^ 1. We additionally require that 
all of the Vjk(t) have the same time dependence. Here we also take Vj k = V k j to be real. 



A. Degenerate 3-level atom 

The degenerate 3-level atom can be solved analytically. Recalling that the Vj k (t) have a 
common time dependence, we choose V(t) = V^t) = a^Vu^t) = P^V^t) = ej Vjj(t). A 
relatively simple solution, presented next, is found by taking €j = 0. More general solutions 
for €j ^ are discussed in an appendix. With ej = and n = 3, Eq(j2J) becomes, 

id x (t) = aV(t)a 2 (t) + /3V(t)a 3 (t) , (3) 
ia 2 (t) = aV{t)ai{t) + V(t)a 3 (t) , 
ia 3 (t) = /?y(t)ai(t)+V(t)a2(t) • 

Now consider the linear combination c(t) = ai(t) + x 0,2 (t) + y % (t), where x and y are 
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some time-independent coefficients that we next determine. Then, 



ic(t) = iai(t) + x ia 2 (t) + y ia 3 (t) (4) 
= aV(t) a 2 (t) + (5V(t) a 3 (t) + x [aV(t) ai (t) + V(t) a 3 (t)} 
+y \J3V(t) ai (t) + V(t) a 2 (t)} . 

Set z = ax + /3y. Then, ic(t) = zV(t)[ ai (t) + (-^L-)a 2 (t) + (-£g-)a 3 (f)]. We require that 
ic(t) = zV(t)c(t). This holds if and only if x — (a + y) / (ax + j3y) and y = ((3 + x)/(ax + /3y). 
After some algebra this leads to the useful cubic equation, 



((3 2 - a 2 )y s + a(2 -a 2 - (3 2 )y 2 + (2a 2 - (3 2 - l)y + a((3 2 - 1) = . 



(5) 



This determines three sets of eigenvalues, {xj}, {yj} and {zj}, and three eigenf unctions, 
Cj (t) = e ~ iz ^\ where A(t) = f*V(t')dt'. Specifically, Cj = E^i-H^* where > 



M 



This matrix, A4, may be inverted, namely, 

/ 



f 1 x x yi ^ 

1 x 2 y 2 
\ 1 x 3 y 3 J 



A 



x 2 y 3 - x 3 y 2 x 3 y x - x x y 3 x x y 2 - x 2 y x 

2/2 - 2/3 2/3 - 2/i 2/i - 2/2 
x 3 — x 2 Xi — x 3 x 2 - Xi 



\ 



(6) 



(7) 



where A = det(M) = x x y 2 + x 2 y 3 + x 3 y 1 - x x y 3 - x 2 y x - x 3 y 2 . 

Now one may express the unperturbed state amplitudes, in terms of the dressed- 

state amplitudes, Cj(t). Using Cj(t) = e - l2 J A (') ; one has a^t) = Y^i=i M-~ l Cj(t) = 
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Ei=iMi le-i *' A(t) - This leads to > 

Pi(t) = K(t)| 2 = -^[{x 2 y 3 - x 3 y 2 ) 2 + (x 3 yi - xiy 3 f + {xxy 2 - x 2 yi) 2 
+2(x 2 y 3 - x 3 y 2 ) (x 3 yi - x x y 3 ) cos ( (21 - z 2 )A(t)) 
+2(x 2 y 3 - x 3 y 2 )(xiy 2 - x 2 y x ) cos({z x - z 3 )A(t)) 
+2(x 3 y 1 - x 1 y 3 )(xiy 2 - x 2 y x ) cos((2: 2 - z 3 )A(t))] , 

p 2 (t) = |a 2 (t)| 2 = -^[(y 2 - y 3 f + (y 3 - Vl f + ( Vl - y 2 f 

+ 2(2/2 - 2/3) (2/3 - C0s((2i - Z2)A(t)) 

+2(2/2 - 2/3) (2/1 - 2/2) cos((^i - zs)4(t)) 
+2(2/3 - 2/1) (2/1 - 2/2) cos((2; 2 - z 3 )A(t))\ , 

•Pa(<) = M^P = [(^3 ~~ X 2) 2 + (Xi - X 3 ) 2 + (x 2 - Xi) 2 

+2(s 3 - x 2 )(xi - x 3 ) cos((zi - «2)A(t)) 
+2(x 3 - x 2 )(a; 2 - x x ) cos((zi - ^)^(*)) 

+2(xi - x 3 )(x 2 - Xl ) cos((z 2 - z 3 )A{t))\ . (8) 

Eq(jHJ) generally gives the transition probabilities to all three final states for arbitrary Vjk(t), 
and corresponds to the more general expression j^J, P&(t) = J2i J2j J ^ki-^kj cos l( z i ~ 
Zj)A(t)}. Since the Xj's, y^s and z/s vary with the Vjk(t), one may then seek conditions on 
the matrix elements Vjk such that the electron populations Pk{to) take any desired values 
at any specified time, t = t . 

B. Complete population transfer 

Next we seek conditions on the V^(i) such that at time t the electron population is 
completely transferred from its initial state i — 1 to a different final state, % = 2 or 3. As 
shown from Eq (jA7|) in an appendix, the condition that 6j = yields a = 2/y — y and 
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(3 = ±1. From the appendix taking (3 = +1 in case i) and using r = ±y one has, 

{x 1: x 2 ,x 3 } = {1,1,-1} , (9) 
{2/i, 2/2, 2/3} = {y+,y-,0} = -{-a + Va 2 + 8, -a - Va 2 + 8, 0} 
= r {m, -n 2 , 0} , 



{zx,z 2 ,z 3 } = {a + y + ,a + —a} = -{a + Va 2 + 8, a - vV + 8, -2a} 
= {-V-, -V+, = r {n 2 , -m, n x - n 2 } , 



A = 2(y_ - y + ) = -2V« 2 + 8 = -2r(ni + ^2) • 

Here ni and n 2 are integers whose values are specified below. This leads to certain allowed 
values of the action integral A(t ) and the relative interaction strength, a, namely, 



3r A(t ) = ±3 J A(t ) = vr , 

V n\n 2 

a = V 12 (t)/V 23 (t) =r (n 2 - ni ) , 

P = V 13 (t)/V 23 (t) = ±1 . (10) 

The same allowed values of A(t ), a and j3 are found in all cases, shown in the appendix. 

One may use the values of the Xj, yj, Zj and A from EqQ in Eq(jSJ) to obtain explicit 
expressions for the transition probabilities, namely, 

Pi(t) = ; - i {n 2 + ^ + n 1 n 2 [l + cos((n 1 + n 2 )rA(t))] 

+ (n 1 + n 2 )[ni cos((2ni - n 2 )rA(t)) + n 2 cos((2n 2 - rii)rA(t)))} , 
P 2W = 777 ! ^{n 2 + ^ + n 1 n 2 [l + cos((n 1 + n 2 )rA(t))] 

— (rii + n 2 )[nx cos((2n! - n 2 )rA(t)) + n 2 cos((2n 2 - ni)rA(t))]} , 

2n x n 2 . 2l ,n 1 + n 2 

P 3W = 7 — T — v sm ( ) rA ^ ' ( U ) 

[nx + n 2 y 2 

where n x , n 2 are odd integers. The constraints on the values of n\ and n 2 are as follows. 
The condition P\ = requires that |(2ni — n 2 ) and |(2n 2 — n\) are both odd integers (i.e. 
tiq and n' D ), and also requires that |(n x + n 2 ) is an even integer. The condition P 2 = 1 
imposes the same requirements. The condition P3 = only requires that 3(711 + n 2 ) be an 
even integer, but P3 > requires that nin 2 > 0. 

Finally a slightly neater solution may now be obtained from the requirements above that 
2rii — n 2 = 3n Q and 2n 2 — n\ = 3n' , where n Q and n' Q are both odd integers (either positive 



or negative). It is easily shown that n\ + n 2 = 3(n + n' Q ) and ri\n 2 = 2n 2 Q + 5n n' + 2n 2 Q . 
Then r = ±(n„ + \n n' + n' 2 )~^, a = r(n' — n Q ) and 3rv4(t ) = 7T. Now, 

= ; ^{n? + + nm 2 [l + cos( (n + n'Jvr (A(t)/A(* )) )] 

+(m +n 2 )[n lC os( n o 7r(A(*)/A(t )) ) +n 2 cos( n' ir (A(t)/A(t Q )) )]} , 

p 2(i) = 777 7 xl^i +n l + ^[1 + cos ( ( n ° + (^(*)M(*o)) )] 

-(n 1 + n 2 )[n 1 cos(n o TT(A(t)/A(t )) ) + n 2 cos( n>(A(t)/A(t )) )]} , 

W) = 7 2 T 2 ^ sin2 ( kno + nXilW/^o)) ) • (12) 
(ni + n 2 J 2 

Here n Q and are arbitrary odd integers, n\ = 2n + n' , and ra 2 = n + 2ra' . At t = to> 
A(t)/A(to) = 1, so that P 2 (to) = 1 with Pi (to) = ^3(^0) = 0. This yields complete population 
transfer from level 1 to level 2 at t — to- We regard n a and n' as the more fundamental 
numbers since they obey the simplest rules. 

III. RESULTS 

Some allowed values of the action integral, A(t ), and the relative interaction strength, 
a, are given in table I. This table includes values of {ni,n 2 }, {n ,n' } and all three {k,k'} 
cases defined in the appendix. From more complete numerical output we confirm that n\ 
and n 2 each acquire all possible odd integer values, although values of the product n\ -n 2 are 
restricted. This is consistent with the condition that m = 2n +n' , and n 2 = n + 2n' , where 
n and n' are arbitrary odd integers. It is also evident that the even integer, n e = n Q + n' , 
takes on all even values. One may also show algebraically that for each value of the even 
integer k in case i) there are two odd values of an odd integer k in case ii), and vice versa. 
The sets of integers {ni,n 2 }, {n ,n' } and {k,k'} are redundant. The three sets of {k,k'} 
correspond to a single set {ni, n 2 }, while the {n , n' } are in one to one correspondence with 
the {n 1 ,n 2 }. 

Numerical calculations for the time dependence of the populations in a degenerate 3-state 
atom perturbed by external interactions with V(t) = V cos(ut) are presented in figures 1- 
4. These results were obtained by using a standard fourth order Runge-Kutta numerical 
integration of EqfJHI) with 6j = and (3 = 1 for various values of a and A(to). For most values 
of a the population transfer into either level 2 or level 3 is always incomplete. A typical 
case is shown in figure 1 where a = 2 and A(t ) = 1.5. We note that P\(t) does return to 
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n\ ■ n 2 


ni 


ni 


n e 


n Q 




n' 


< 


n e 




a 








k 


-k' 


k 


k' 


-k 


k' 






5 


±1 


±5 


±2 


±1 


±1 


±3 


±3 


±2 


±1.656 


T2.530 


5 


±5 


±1 


±2 


±3 


±3 


±1 


±1 


±2 


±1.656 


±2.530 


9 


±3 


±3 


±2 


±1 


±1 


±1 


±1 


±2 


±2.221 


0.000 


11 


±1 


±11 


±4 


±3 


±3 


±7 


±7 


±4 


±2.456 


T4.264 


11 


±11 


±1 


±4 


±7 


±7 


±3 


±3 


±4 


±2.456 


±4.264 


17 


±1 


±17 


±6 


±5 


±5 


±11 


±11 


±6 


±3.053 


T5.488 


17 


±17 


±1 


±6 


±11 


±11 


±5 


±5 


±6 


±3.053 


±5.488 


23 


±1 


±23 


±8 


±7 


±7 


±15 


±15 


±8 


±3.551 


T6.487 


23 


±23 


±1 


±8 


±15 


±15 


±7 


±7 


±8 


±3.551 


±6.487 


27 


±3 


±9 


±4 


±1 


±1 


±5 


±5 


±4 


±3.848 


Tl-633 


27 


±9 


±3 


±4 


±5 


±5 


±1 


±1 


±4 


±3.848 


±1.633 


29 


±1 


±29 


±10 


±9 


±9 


±19 


±19 ±10 


±3.988 


T7.353 


29 


±29 


±1 


±10 ±19 


±19 


±9 


±9 


±10 


±3.988 


±7.353 


35 


±1 


±35 


±12 Til 


±11 


±23 


±23 ±12 


±4.381 


T8.128 


35 


±5 


±7 


±4 


±1 


±1 


±3 


±3 


±4 


±4.381 


T0.478 


35 


±7 


±5 


±4 


±3 


±3 


±1 


±1 


±4 


±4.381 


±0.478 


35 


±35 


±1 


±12 ±23 


±23 


±11 


Til ±12 


±4.381 


±8.128 



TABLE I: Some values for the action integral, A(to), and the relative interaction strength, a, 
allowed for total population transfer. These values are found using Ea (|10j) subject to the conditions 
listed. Here n e = n Q ± n' Q , n Q = ^(2ni — 712) and n' Q = \{2ni — rii). 

zero twice in each period, T, of the external field, but that complete population transfer to 
an initially unoccupied state never occurs. We also ran calculations (not shown) for a = 8.5 
and A(t ) = 5.8, which are not values that lead to complete population transfer from Eq()10|). 
Again complete transfer to an initially unoccupied level never occured. However, there were 
rapid oscillations in the populations of all states except near t = T/4 and t = 3T/4, where 
none of the populations oscillated rapidly. This was similar to figure 4. This appears to 
correspond to the onset of complete population transfer, which occurs at t = T/4 and odd 
multiples of T/4, as seen in figures 2-4. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Time (in units of the field period, T=2TC/&)) 



FIG. 1: Occupation probabilities as a function of time over one period, T, of the external field, 
V(t) = Vocos(a;i). The solid line is the probability, P%(t), that the electron is transferred to 
the target state, namely state 2. The long dash line is the probability, P\(t), that the electron 
is in state 1, the state in which it began initially. The short dash line is the probability, Ps(t), 
that the electron is in state 3. In this figure a = V12/V23 = 2 and the action area at t = to 
is A(to) = J *° V(t')dt' = 1.5. Since a and A(to) do not correspond to values given by Ea (|lU|) . 
complete transfer to level 2 does not occur. 




Time (in units of the field period, T=2TC/(B) 



FIG. 2: Occupation probabilities as a function of time. The solid line denotes i^(£); the long 
dash line denotes Pi(t); and the short dash line denotes Ps(t). In this figure we use the allowed 
values, a = P12/V23 = and A(to) = j Q ° V(t')dt' = 2.221, corresponding to n\ = 3 and ri2 = 3 
(n = n' = 1) in Eq^nj). Complete transfer to level 2 from level 1 occurs at t = to = T/4 and 
again at odd multiples of to. 



9 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Time (in units of the field period, T=2TC/(B) 



FIG. 3: Occupation probabilities as a function of time. The solid line denotes P2(t); the long dash 
line denotes Pi(t); and short dash line denotes P3(t). In this figure a = —2.530 and A(to) = 1.656, 
corresponding to n\ = 1 and ti-2 = 5 (n = — 1, n' = 3) in Eq^UJ). 

Calculations using a few values of a that permit complete transfer to level 2 are shown in 
figures 2-4. Our numerical codes give the same results as the analytic expressions of Eq ffTTj) . 
This provides a check that our algebra is correct. We note that 2ni " 2 , 2 < 1 and that the 
maximum value of Ps(t) occurs for m = ri2, where P3 max = |, consistent with figure 2. This 
corresponds to a = so that direct transitions from level 1 to level 2 are forbidden. Transfer 
to level 2 occurs via the intermediate level 3. Transfer from level 1 to level 2 and back is 
complete, and occurs periodically. In general a = occurs when n\ = n 2 = 3n odd , where 
n dd is any odd integer. This appears to give the simplest condition that allows complete 
population transfer. In this case the action area is A(t ) = n oc id^/V2- 

Calculations for two other values of a that allow complete transfer to level 2 are shown in 
figures 3 and 4. We see that complete transfer occurs twice in one period of the oscillating 
field but that the frequency of the "side bands" increases as aA increases. We note that 



a = ±w^- (n 2 — ni) becomes either large (n 2 ni, or vice versa), or small (ni ~ n 2 ) 



as n\n 2 increases, while A(to) increases as \Jn 2 n 2 /2. When complete population transfer 
occurs, the population lingers in level 2, as seen in figures 2-4. It can be shown (2^ that 
near t = t 1 — P 2 (t) varies as [u(t — to)} 4 . 

Additional control 25] may be achieved by changing the shape of V(t). This can be 
used to control how long the population remains near unity in level 2, for example. An 
interesting example is the case of an ideal, sudden 'kick' produced by V(t) = A 5(t — 1 ). If 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Time (in units of the field period, T=2TC/(B) 



FIG. 4: Occupation probabilities as a function of time. The solid line shows P 2 (t); the long dash line 
shows Pi(t); and the short dash line shows Ps(t). Here a = 8.128 and A(to) = 4.381, corresponding 
to m = 35 and 712 = 1 (n = 23, n' = —11) in Ea lfTU)l . Since these values are allowed, the electron 
population is completely transferred to level 2 from level 1. 

^4o ~~ > n-odd^/V^ for example, then the mere presence of level 3 allows transfer from level 1 
to level 2 without any direct transfer from level 1 to level 2. With this ideal 'kick', level 2 
is unoccupied before t = t , and fully occupied after t = t . 

To obtain analytic solutions to Eq.(J2J) for a 3-level atom, we imposed degeneracy on 
the 3-level manifold, e.g. by taking Ej = E = 0. If E ^ 0, then E can be removed 
from Eq ([T|) by an overall phase transformation, e lEt . When the levels are not degenerate 
Uij = (Ei — Ej)/h 0, and the population transfer is incomplete. The effect of this 
population leakage can be calculated by expanding the transition amplitude 02 in a power 
series in time with terms small in lo^t. After some algebra we obtain from EqQ the 
difference between the degenerate and nearly degenerate populations of level 2, namely, 
AP 2 (t) ~ t^[2(2cj 13 - cj 12 )l / i2(0)\/i3(0)y23(0) + ul 2 V? 2 (0)]t 4 , where 1^(0) denotes the value 
of Vij(t) at t — 0. By the time the occupation probability P 2 reaches its first maximum at 
t Q = T/4 = using 0-^12(0) = p-'VniO) = V 23 (0) = V(0) = ±^f^u, the difference 
becomes, AP 2 (t ) ~ ^f[^ ni n 2 (n 2 - n 1 )(2^ - ±f) + (n 2 - m) 2 (^) 2 ]. The choice 
of the field frequency, u, involves a trade-off between the duration of time the population 
remains in state 2 and the population leakage. The higher the frequency, u, the smaller 
the population leakage, AP 2 , but the shorter the duration time, T s , that the population 
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remains in level 2. This effect of population leakage is similar as that for a 2-level atom J^. 
Numerical calculations indicate that this population leakage grows rapidly in time. 

IV. DISCUSSION 

In this paper we have considered transitions amplitudes and probabilities in a 3-level atom 
with degenerate states. We also assumed that the interaction matrix elements, Vij(t) all have 
a common time dependence. In realistic atomic systems, energy levels are seldom, if ever, 
exactly degenerate. Also levels outside 3-level manifold usually exist. Consequently, use of 
our results is restricted to external fields with frequencies in the range, io m % n < uj < u max . 
Here hj min is the energy splitting of the nearly degenerate states, and Hu> max is the energy 
difference between the 3-level manifold and the closest state in energy outside the manifold. 

It is instructive to compare our 3-level results to simpler 2-level results Q Q The 
equations for a 2-level atom may be recovered from our 3-level equations by taking Vyi and 
V13 to zero, which corresponds to a — > oo (with an y \ej /a\ < oo), as may be understood 
from Eq(j3J. For the 2-level atom it has been shown 25|,|26( that aj(t) may be expressed in 

2 | 1 y + 

terms of eigenstates Cj(t) via the relation, Cj(i) = z2j=i MijijWi where M.^ = I 

^ ' \ 1 y- 

Here y± are the same as those given by Eq (jA8|) . This matrix may be inverted to give 

the population amplitudes, a\{t) = i[y_ci(t) — y^c^it)] and a 2 (t) = -A— c\{t) + c 2 (t)], 



where A = 2y 1 + ( e2 2 £1 ) 2 is the determinant of M.ij. A simple solution occurs when 
e x = e 2 = e, namely, a x {t) = e~ ieA{t) \[e~ iA ^ + e +iA W] = e~ ieA ^ cos (A(t)) and a 2 (t) = 
e -ieA(t) i[ e -iA(*)_ e +iA(t)] = _ ie -ieA(t) s i n (A(t)). Note that Pi = 1-P 2 . One may then simply 
determine the conditions under which the population of level 2, P2(to), takes on any desired 
value at time t . If e± ^ e 2 , it is easily shown that i^o) < 1/(1 + (e 2 — ei) 2 /4) < 1. Thus the 
diagonal matrix elements of Vij(t) prevent complete transfer to the initially unoccupied level 
when ei and e 2 are unequal, and, when they are equal, they simply contribute an overall 
phase. When ei = e 2 any value of P 2 (to) between and 1 can be found. In particular, 
-^2(^0) = 1 if A(to) = Jq° Vi2(t')dt' = noddTt/Z. This allowed value of A(t Q ) for complete 
transfer to level 2 differs from the allowed values given in Eq ljlOj) for degenerate 3-level 
atoms. By comparison to a 2-level atom the conditions for complete transfer are generally 
more complex in a 3-level atom, as discussed in the appendix. For example, complete 
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population transfer can occur in a 3-level atom when e 2 7^ £3. In both cases, however, 
the duration of time spent in the transferred state can be controlled by adjusting the time 
dependence of the external field. Also if the states are not quite degenerate then population 
leakage occurs. For a harmonic V(t) the population leakage in a 2-level atom varies as 
A-Bjf^o) — ^(^) 6 (ui 2 /uj) 2 . For 2s — 2p transitions in atomic hydrogen complete population 
control can be nearly achieved in this manner using a radiation field with wavelengths (and 
intensities) ranging from a few fim (with about 10 12 W/cm 2 ) up to a few cm (with about 
10 4 W/cm 2 ). At wavelengths below a few fim coupling to nearby n = 3 atomic levels can 
be significant, and this 2-level model breaks down. 

While we have not provided calculations in this paper for specific experiments, some 
general guidelines for experimental tests and applications are evident. First, the model must 
be valid, so that the frequency of the external interaction is limited by io m m < uj < u max . 
If the external interaction varies harmonically, V(t) = Vocos(c<jt), then the first line of 
Eq ljlO)) imposes a constraint between Vq and uj, namely, Vq/oj = ±*/^p^. For a harmonic 
interaction the duration of the time, T s , that the population remains in level 2 varies inversely 
with uj; the higher the frequency, uj, the smaller the time the population remains in level 2. 
Our degenerate 3-level model can be applied to systems with dipole selection rules. A dipole 
selection rule can correspond to a = 0. Transfer from level 1 to level 2 (which is the dipole 
forbidden transition) is complete at time to if the two allowed transition matrix elements, 
V23 and Vi 3 , are equal in absolute magnitude and A(t ) = n dd'n/2. In some atomic systems 
state 2 could decay, with a lifetime, Td, to another level outside the degenerate manifold and 
be lost. Such loss can be controlled by adjusting uj. Our model also rests on degeneracy. If 
the three levels are not degenerate the transfer of population to level 2 is incomplete. We 
call this population leakage. As discussed above this population leakage may be minimized 
by using high frequency external interactions, i.e. uj ^s> uj min , but at a cost to duration of 
time, T s , the population remains in level 2. 

Using nearly degenerate states, the duration time, T s , can be further controlled 
by adjusting the shape of V(t). As mentioned in section III for example, the transfer is 
complete, instantaneous and permanent when V(t) = n od d^S(t — t ). Such a quick, hard 



pulse is called 27|, |28|, [29] a 'kick'. Two practical limitations on this ideal model are the 



impossibility of producing a signal that varies as 5(t — to), and the existence of an infinitely 
wide spectrum of high frequency components with frequencies uj > 0J max in the Fourier 
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spectrum of 5(t — to). Fortunately these two difficulties can both be addressed by using 
'kicks' of finite width in time. In some cases it may be possible to design a 'kick' so that 
its duration is short compared to any other changes in the system, so that a finite 'kick' 
may be sensibly represented by 5(t — t ). If, in addition, the energy levels outside the 
(nearly) degenerate manifold have a large energy gap hjj maxi then it may be possible that 
ujmin < ijj < u max and our model may be applicable. Applying V(t) ~ 5(t — t ) to nearly 
degenerate atomic systems leads to the 'gedanken' question, what happens when one tries to 
force the transition to occur within a small time interval about to in a degenerate quantum 
system where At is large? 

Let us next consider some mathematical features that occur when degenerate states are 
used. The exact probability amplitude for a transition from an initial state \i) to a fi- 
nal state \ f) may be expressed in terms of the time evolution operator U(ti,tf), namely, 
a fi = (f\U(t f ,ti)\i), where (setting U = and t f = t), U(t,0) = Te^ti^** = 
^y- f* ... J t TV(t n )....V(ti)dt n ...dti. Here T is the Dyson time ordering operator which 
places the interactions in the order of increasing time so that for t > t', TV(t)V(t') and 
TV(t)V(t') = 0, and vice versa for t! > t, whence T[V(t),V(t')} ^ 0. In any intermediate 
time step from tj to t j+1 , the matrix element is J * e ( - El ~ El '^ < I\V(tj)\I' > dtj. Without 
the phase e < ^ El ~ El '" >tj there is no time ordering since T[V(t), V(t')] — > as T — > 1. If all the 
states are degenerate, then AE = Ej — Ef — » in every time step. This corresponds to 
taking T -> 1 in U(t, 0) = Te'^o v ( t ') dt ' . In this limit of degeneracy, U(t, 0) e~ l ti v ( t ') dt \ 
corresponding to the approximation of Magnus j^, where [V(t), V(t')] — > 0. In this limit 
there are no time correlations in the time propagation. Imposing degeneracy removes time 
ordering. Moreover in standard stationary scattering with H = H + V, the Fourier trans- 
form of e l fo( E - H "> dt is irc5(E — H), which is the energy conserving part of the Green function, 
G = E-H±ir] = iwfiiE ~ H) =F Pv^zjj- The energy non-conserving quantum fluctuations, 
P v E ^_ H , imposed by the asymptotic ±irj initial condition, carry the direction of time, i.e., 
outgoing or incoming scattered waves. Here P v -^-^ is the principal value integral that omits 
the contribution at E = H. The Fourier transform of P v E i _ H is sign(t) = T — 1. Degener- 
acy, P v F ^_ H — > 0, again removes time ordering. Thus, imposing degeneracy mathematically 
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corresponds to, 



AE -f <=^> [ e {El ~ E ^ < I\V(tj)\r > dtj -> f < I\V(tj)\I' > dtj 
Jo Jo 

^ T[V (t),V(lf)] -> £/(t,0) -> e^^ v ^ 



e i(E-H )t Te -if*V(t')dt' _^ e iJ*(E-H)dt 
1 1 P 1 



E-H±irj E-H E-H 
T -> 1 . (13) 

As a consequence of using degenerate states, quantum energy fluctuations in intermediate 
states are eliminated. The Hilbert space is restricted so that At = H/AE — > oo. In addition 
the minimal size of the wavepacket, A£ = vAt, now extends to infinity, so that the quantum 
system may not be decoherently decoupled from its macroscopic environment. This effect is 
evident in 2s — 2p transitions in hydrogen caused by the impact of ions, where the range of 
the interaction goes to infinity as the levels become degenerate Within the interval 

At time cannot be localized, e.g. reproducibly observed. Time ordering i.e. causality, is 
not operative and there is no 'flow of time'. In the classical limit as At — > 0, time can be 
measurably localized, being defined by two closely spaced decoherent observable events. In 
this classical limit causality strictly holds. 

We note that at sufficiently high uj all bound levels in the atomic elements become nearly 
degenerate. Hence, if one can deal with high energy continuum states, our approach might 
be useful in applications involving fourth generation synchrotrons that produce intense high 
frequency fields. An extension of this method might also used for control transitions in high 



M, I35J- RWA is used to describe 



Rydberg states |33j, including adiabatic rapid passage 
coherent storage of information in photonic states The approach developed here may 
be useful for modeling information transmission and storage in atomic states Q| in a new 
way. 



V. SUMMARY 



We have shown that in a 3-level atom with degenerate energies, electron population is 
completely transferred via an external interaction at a designated time, to, from an ini- 
tially occupied level (level 1) to a designated initially unoccupied level (level 2) under two 
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conditions. The first condition for complete transfer is that the ratio of the matrix ele- 
ments of the external interaction, V^i), satisfy V\i(t) /Vzzif) = a = ±^J ' -^^{n 2 — Tlx), and 
Vi3(t)/V 23 (t) = (3 = ±1. The second condition is that at t — to the action area of V(t) 
satisfy A(t Q ) = jj° V(t')dt' = iy/^pf , where we have set V 23 (t) = V(t). Here m and n 2 
are integers such that n\ = 2n + n' and n 2 = n + 2n' , where n Q and n' are any odd integers. 
The duration of time the transferred population remains in level 2 can be controlled either 
by varying the frequency, u, of the external potential, or by varying the shape of V(t). 
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APPENDIX A: METHOD TO FIND CONDITIONS FOR COMPLETE POPU- 
LATION TRANSFER 

Here we show how to find the conditions on the Vjk such that the electron population is 
completely transferred from its initial state i — 1 to a different final state, % = 2 or 3. A 
convenient way to begin is to take dPk(t)/dt = at t = to- This gives both maxima and 
minima for P*.. Using Eq(|5|) this yields, 

sin((^ - Zj)A(t )) = (Al) 
(zi - Zj)A(t ) = m i: j-K (rriij ^ 0) (i, j = 1,2,3) . 

Using cos((zj — Zj)A(t )) = cos(my7r) = (— l) miJ , one has, 

Pi = |«l| 2 = ^[(^22/3 ~~ X 3V2) 2 + ( X 3Vl ~ x lll3) 2 + ( x lV2 ~ x 2Vl) 2 

+2(-l) mi2 (x 2 y 3 - x 3 y 2 )(x 3 y 1 - x x y 3 ) 
+2(-l) mi3 (x 2 y 3 - x 3 y 2 ){xiy 2 - x 2 y x ) 
+2(-l) m2a (x 3 y 1 - xiy 3 ){xiy 2 - x 2 yi)\ 
= ^{{-l) k ^x 2 y 3 - x 3 y 2 ) + (-l) fc2 (x 3 yi - x iy3 ) + (-l) ks (x iy2 - x 2Vl )} 2 , (A2) 
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where k± + k 2 = mi 2 , k\ + k 3 = mi 3 and k 2 + k 3 = m 23 . Now, 

k = h = \{m X2 + m 13 - m 23 ) = \{ ^~ (zi - z 2 ) + ^-^-(zi - z 3 ) - ^ L ^-(z 2 - z 3 )) 

2 2 71 71 71 

= ^(*-*) 

7T 

&2 = ^(»"12 + m 23 - "lis) = J(— —( z l - z 2) + ^~^-(z 2 ~ Z 3 ) - ^^-(Zi - Z 3 )) = 
2 2 IT 71 71 

k' = k 3 = -(mi3 + m 23 - m X2 ) = -( — — (21 - z 3 ) H — (z 2 - 23) —{zx - z 2 )) 

2 2 71 7T 7T 

= ^)(^-*3) (A3) 

71 

This yields extrema for P^ at t — to, namely, 

Pi = ^ [ (-1)^22/3 - x 3 y 2 ) + (x 3 yi - x x y 3 ) + (-l) fc '(xi?/ 2 - x 2 ?/i)] 2 

^2 = ^[(-l) fe (2/2-y3) + (y3-yi) + (-l) fc '(yi-Z/2)] 2 

P 3 = -^[(-l) k (x 3 - x 2 ) + (xi - x 3 ) + (-l) fc '(a; 2 - x,)] 2 (A4) 



where and k! are non zero integers that are not both even, 
a. Conditions for complete transfer 

We now seek conditions such that at t — to, P\ — 0, P 2 — 1, and P 3 = 0. There are three 
sets, k'}, that satisfy these conditions, namely, i) k is even and k! is odd, ii) k and are 
both odd, and iii) k is odd and k' is even. 

Case v . k is even and k! is odd. Then, 

p i = T^K 2 ^ - ^2) + {x 3 y x - x x y 3 ) - {x x y 2 - x 2 y x )f 



A 2 

1 4 

p 2 = ^2 [(^ 2 ~~ y*) + ~ ~ fa " y^ 2 = ~ y^ 2 
1 4 

^3 = ^[(^3 - x 2 ) + (xi - x 3 ) - (x 2 - x^] 2 = ^(^1 - x 2 ) 2 . (A5) 

From the condition that P 2 = 1 one has [2(y 2 — y^/A] 2 = 1, where ?/! 7^ y 2 , so that 
A = ±2(y 2 — The condition that P 3 = gives [2(xi — x 2 )/A] 2 = 0, which yields 
X\ — X 2 — X. Using this result one then obtains P\ = [(x + x 3 )(yi — y 2 )/A] . Requiring 
Pi = 0, one has x 3 = —x. Now, using x x = x 2 = —x 3 = x, we have from below Eq(J7|), 
A = 2x(y 2 — yi). Finally since A = ±2(y 2 — yi), one has x\ = x 2 = —x 3 = x = ±1. 
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Case ii : k and k' are both odd. Then, 

Pi = -^[-(^3 - %m) + {xzVi - - (%m - x 2 yi)] 2 

P2 = M-(v2 - a.) + (vs - vi) - (vi - 1/ 2 )] 2 = Mv*~ ( A6 ) 

■P3 = a*[-0&3 - ar 2 ) + (xi - ac 3 ) - (x 2 - x x )f = ^(x 1 - x 3 ) 2 . 

From the condition that P 2 = 1 one has [2(7/3 — 2/i)/A] 2 = 1, where j/i 7^ 1/3, so that 
A = ±2(2/3 — yi). The condition that P3 = gives [2(xi — Xs)/A} 2 = 0, which yields 
x\ = x 3 = x. Then one obtains Pi = [(x + x 2 )(yi — 1/3)/ A] 2 . When Pi — one has 
x 2 = —x. Now, using xi = —x 2 = x 3 = x, we get from below Eq©, A = 2x(yi — 2/3). Since 
A = ±2(y 3 — ?/i), one has Xi = —x 2 = x 3 = x = =Fl. 

Case in : k is odd and k' is even. This is similar to case i). One obtains, A = 2x(y% — y 2 ) 
and —x\ — x 2 — X3 = x = ±1. 

Complete population transfer, i.e. P 2 — 1, requires that x = ±1, so that the relative 
interaction strengths a and (3 must satisfy, 

2 e 2 - ei , . 

a = 2/ ± 5— ± (£3 - e 2 ) 

2/ 2/ 

= ±1 + f?JLf* (A7) 
1/ 

Next we use these results to find the conditions on V^- required for complete population 
transfer to occur. 

b. Quantum numbers for A(to) and a 

We now note that relatively simple conditions occur when ti = e 2 = e 3 = e, namely, 
a = ^ — y and (3 = ±1. In this case Eq (|A7|) reduces from a cubic to an easier quadratic 
equation in y. This quadratic equation corresponds to three roots of the more general cubic 
equation, namely, 

fa} = fr + ,y-o} = r n + f* TI .^g^.o} . (A8) 

Now for case i) when k is even and k 1 is odd, one has from above, x = {±1, ±1, =Fl} 
and y = {y + ,y_,0}. From Eq (|A8|) . z = {=FJ/-, Ty+, T®}- Then, with E = A(to)/n, the 
useful integers k and k' are: k = ^&l(zi — z 2 ) = ±E\/ a 2 + 8, and k! = ^^-(z 2 — Z3) = 
±f (3a - \Ja 2 + 8). After a little algebra one obtains (k - k')(2k + k') = 18E 2 = n x n 2 . 
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For case ii) when k and k' are both odd, x = {±1,^=1, ±1} and y = {y + ,0,y~}. and 
z = { Ty _, T a,Ty+}- Then, k = «&( Zl - z 2 ) = ±f (3a + and k' = ^(z 2 - 

z 3 ) = (3a - Va 2 + 8). One may now verify that (2k + k')(k + 2k') = 18E 2 = n x n 2 . 
For case iii) when k is odd and k! is even, x = {=Fl,±l,±l} and y = {0,y + ,y_}, and 
z = {T«, TV-, TV+}- Then, k = «&( Zl -z 2 ) = T^a+V^TS), and k' = *M(z 2 -z 3 ) = 
±E\/ a 2 + 8). One may again verify that (2k 1 + k)(k' — k) = 18E 2 = ri\n 2 . 

Using a 2 = k 2 /E 2 — 8 one readily finds in case i), where k is even and k' is odd, that 
k = |(ni + n 2 ) and k! = \(n 2 — 2ni). In case ii), where and k! are both odd, one has 
k = g(2ni — ^2) and fc' = \(2n 2 — n\). In case iii) k = |(m — 2n 2 ) and fc' = |(rii + n 2 ). In 
all three cases one readily obtains the same equation for the allowed values of A(t ) and a, 
namely, 

a = ±J-^— (n 2 -n x ) . (A9) 
V n x n 2 

It is remarkable that all three cases give the same formula for complete population transfer. 
The top equation corresponds to a constraint on the value of the action, J Q * V(t')dt', i.e. the 
"area" of the external interaction, V(t), or the phase shift caused by the external interaction. 
The bottom equation above gives the relative values of the V^-. 



APPENDIX B: GENERAL SOLUTIONS FOR THE DEGENERATE 3-LEVEL 
ATOM 

In the main text a solution for the degenerate 3-level atom was found in the case when 
tj = 0. We now choose e^Vj^t) = V(t) = V 23 (t) = a~Vi 2 (0 = (3~ l V n (t). Then Eq0 
becomes, 

ia^t) = eiV(t)ai(t) + aV(t)a 2 (t) + (3V(t)a 3 (t) 
ia 2 (t) = aV(t) ai (t) + e 2 V(t)a 2 (t) + V(t)a 3 (t) 

ia 3 (t) = pV(t)ai(t) + V(t)a 2 (t) + e 3 V(t)a 3 (t) . (Bl) 
Consider a linear superposition c(t) = a x (t) + x a 2 (t) + y a 3 (t), where x and y are some 
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time-independent coefficients that we next determine. Then, 

ic(t) = iai(t) + x ia 2 (t) + y ia 3 (t) 

= dVit) ai {t) + aV(t) a 2 (t) + (3V{t) a 3 {t) 
+x[aV(t) ai(t) + e 2 V{t) a 2 (t) + V(t) a 3 (t)} 
+y[(3V{t) ai (t) + V(t) a 2 (t) + e 3 V(t) a 3 (t)} . (B2) 

Set z = e l + ax + (3y. Then, ic{t) = zV(t)\ ai (t) + (5^)a 2 (t) + (|ggjL)a 3 (t)]. 
We require that ic(t) = zV{t)c{t). This holds if and only if z = e\ + ax + x = 
(a + e 2 x + y) / (ei + ax + (3y) and y = {(3 + x + e 3 y) / (ei + ax + . After some algebra this 
leads to the useful cubic equation, 

[((3 2 - a 2 ) + a(3(e 2 - e 3 )]y 3 + [a(2 -a 2 - (3 2 ) + (3(26, -e 2 - e 3 ) + a(e x - e 3 )(e 2 - e 3 )]y 2 
+ [(2a 2 -{3 2 -l) + a(3(2e 3 -e x - e 2 ) + (e x - e 2 )( ei - e 3 )\y + [a(/3 2 - 1) - /3( ei - e 2 )\ 
= . (B3) 

This determines three sets of eigenvalues, {xj}, {yj} and {zj}, and three eigenf unctions, 
Cj (t) = e - i2 ^W, where A(t) = f*V(t')dt'. 

When 6j 7^ complete population transfer can also occur. Using ej = e ^ adds 
nothing new, since it produces only an overall phase e let . A slightly more general, but 
less elegant, solution for the degenerate 3-level equations is found when ei = e 2 ^ e 3 . 
Again in this case a = | — y and (3 — ±1. The quadratic equation for y now becomes 

2 a(1 _ Q 2 )+ ( )2±(ei _ e3) _ 2 = ^ = a( l- a 2 )+ ( 3) 2 ±(ei _ e3) 

^ 1— o^±ct(ei— £3) f 1— a^±a(e\— £3) 

same as that one above, except that a — > <x This includes an additional parameter, (ei — e 3 ), 
gives simple allowed values of a, and imposes a constraint between a and (ei — e 3 ). More 
complex conditions for complete population transfer might exist when e\ 7^ e 2 , in which 
case Eq (1B3j) is cubic in y. Finally it might be possible that solutions exist when the matrix 
elements Vij(t) have different time dependences, i.e. when a and (3 depend on t. 

Conditions for complete transfer into level 3 at t — to may be generally found similarly 
by interchange of the indices 2 and 3, with corresponding interchange of a with (3 and x 
with y. 
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